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Evolutionary structure searches are used to predict stable phases with unique stoichiometrics in the hydrogen- 
rich region of the magnesium/hydrogen phase diagram under pressure. MgH4, MgHi2 and MgHi6 are found 
to be thermodynamically stable with respect to decomposition into MgH2 and H2 near 100 GPa, and all lie on 
the convex hull by 200 GPa. MgH4 contains two H~ anions and one H2 molecule per Mg 2+ cation, whereas 
the hydrogenic sublattices of MgHi2 and MgHi6 are composed solely of H^ - molecules. The high-hydrogen 
content stoichiometrics have a large density of states at the Fermi level, and the T c of MgHi 2 at 140 GPa is 
calculated to be nearly three times greater than that of the classic hydride, MgH 2 , at 180 GPa. 
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I. INTRODUCTION 

In 1935, Wigner and Huntington predicted that hydrogen, 
which exists in the paired molecular state at ambient pres- 
sures and temperatures, would become an alkali-metal-like 
monoatomic solid when compressed to pressures exceeding 
25 GPa 0. This turned out to be a bit of an underes- 
timate. Structure searches based upon density functional 
theory calculations have identified a number of molecular 
(P < 400 GPa) Q and quasi-molecular or atomic (P = 
0.5 — 5 TPa) (3][4l structures in the cold phase diagram, and 
experiments in diamond anvil cells show that the insulating 
phase III with paired hydrogens is stable over a broad tem- 
perature range and up to at least 360 GPa Q. Recent ex- 
perimental work at room temperature above 220 GPa, which 
showed a pronounced softening of one of the Raman active 
molecular vibrons upon compression [6], and potential con- 
ductivity [ 7 ] has generated much excitement (8). The newly 
discovered phase IV of hydrogen is thought to be a mixed 
structure composed of layers of molecular units, as well as 
weakly bonded graphene-like sheets 0[lOl. Superconductiv- 
ity at pressures of 450 GPa and temperatures up to 242 K has 
been predicted in the molecular phase [11], whereas the super- 
conducting transition temperature, T c , may approach 764 K 
for monoatomic hydrogen near 2 TPa lfT2ll . 

What are the structural motifs that compressed hydrogen 
adopts when doped with an electropositive element? Theoret- 
ical work has predicted the presence of hydridic H~ atoms, 
H2 - molecules, H~ • • • H2 motifs, symmetric H3 molecules, 
polymeric (H3 ) oc chains, and sodalite cage structures, with 
the nature of the hydrogenic sublattice depending upon the 
identity of the alkali or alkaline earth metal, and the pressure 
fl3l420l . The latter clathrate-like cage which encapsulated 
the calcium cation in CaH 6 was shown to be susceptible to a 
Jahn-Teller distortion, giving rise to a remarkable electron- 
phonon coupling parameter of 2.69 at 150 GPa with a con- 
comitant T c of -225 K 11201 . 
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Herein, evolutionary structure searches are used to seek 
out the stoichiometries and structures of magnesium polyhy- 
drides, MgH n with n > 2, under pressure. MgH4 is pre- 
dicted to become stable with respect to decomposition into 
MgH2 and H2 near 100 GPa, and it remains the most stable 
stoichiometry until at least 200 GPa. Cracra-MgH 4 contains 
two hydridic hydrogens and one hydrogen molecule per Mg 2+ 
cation, and it becomes metallic as a result of pressure-induced 
band overlap. However, due to the low density of states at 
the Fermi level, the estimated T c at 100 GPa is only -10 K 
higher than that of P63 /mmc-MgH 2 at 180 GPa. Two other 
stoichiometries, MgHi2 and MgHi6, are calculated as being 
thermodynamically stable at higher pressures. Because these 
hydrogen rich phases exhibit an "H^ - belt" surrounding the 
Mg 2+ cations, and do not contain any hydridic hydrogens, 
they have a high density of states at the Fermi level. Assum- 
ing typical values of the Coulomb pseudopotential, the T c of 
MgHi2 is calculated as being between 47-60 K at 140 GPa. 
Phases with a greater ratio of H^iH - tend to have a larger 
density of states at the Fermi level, and a higher T c . 

II. COMPUTATIONAL METHODS 

The structural searches were performed using the open- 
source evolutionary algorithm (EA) XtalOpt release 8 along 
with the default parameter set l2T1l . Evolutionary runs were 
carried out on the MgH 2 stoichiometry at 0, 25, 50, 100, 
150, 200 and 250 GPa with cells containing 2, 3, 4 and 8 
formula units (FU). At 200 GPa additional searches on cells 
with 5 and 6 FU were performed. Only the known a, e and 
(-phases were recovered, with the P6s/mmc structure re- 
maining the most stable up to the highest pressures consid- 
ered. Exploratory searches at 200 GPa revealed that MgH n 
with n = 2.5,3,3.5,4.5,5,... were noticeably less stable than 
those with even n. Thus, a more refined search was carried 
out at 100 and 200 GPa which was restricted to even n rang- 
ing from 4-16 with 2-4 FU for MgH 4 and MgH 6 , 2-3 FU for 
MgHg and 2 FU otherwise. Additional searches for n = 4 
and n = 6 were carried out with 8 and 5 FU at 100 and 
200 GPa, respectively. Duplicate structures were detected us- 
ing the XtalComp [ 22 ] algorithm. The spglib package l23l 
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was used to determine space-group symmetries. 

Geometry optimizations and electronic structure calcula- 
tions were performed using using density functional theory 
(DFT) as implemented in the Vienna ab-initio simulation 
package (VASP) (24). The exchange and correlation effects 
were treated using the Perdew-Burke-Ernzerhof (PBE) func- 
tional 1 25 ] with plane-wave basis sets and a kinetic energy 
cutoff of 600 eV. The hydrogen Is 1 and magnesium 2p 6 3s 2 
electrons were treated explicitly and the projector- augmented 
wave (PAW) method (26) was used to treat the core states. The 
/c-point grids were chosen using the T-centered Monkhorst- 
Pack scheme. The number of divisions along each reciprocal 
lattice vector was chosen such that the product of this number 
with the real lattice constant was 40 A for the final geome- 
try optimizations, as well as at least 50 A for the electronic 
densities of states (DOS) and band structures. 

In situations where band gap closure occurs as a result 
of pressure induced broadening, and eventual overlap of the 
valence and conduction bands, standard density functionals 
predict too low metallization pressures. For this reason we 
have calculated the DOS of select structures using the HSE06 
screened hybrid functional, which has been shown to give 
good accuracy for band gaps (27). Due to the immense com- 
putational expense involved in hybrid calculations, the ge- 
ometries employed have been optimized with PBE. Recently, 
it has been shown that inclusion of Hartree-Fock exchange 
in functionals such as HSE06 or PBEO can have a significant 
impact upon the calculated transition pressures between dif- 
ferent phases (28] [29). Moreover, a study considering liquid 
nitrogen at pressures up to 200 GPa (and finite temperatures) 
has shown that the structural relaxations with a hybrid func- 
tional can lead to large Peierls distortions (which inevitably 
have an impact on the band gap) (30). It has been proposed 
that nontrivial exchange correlation effects can become quite 
important at extreme pressures, in particular when electron 
localization occurs ("electride" behavior (3D ), and/or when 
the semicore electrons interact (29) . Since we do not observe 
significant broadening of the Mg 2p bands at the pressures 
considered here, it may be that structural relaxation with a hy- 
brid functional will yield similar results as PBE. This will be 
considered in future studies. However, as expected, our calcu- 
lations show that hybrid functionals increase the pressures at 
which band gap closure is predicted to occur. 

Phonon calculations were performed using VASP com- 
bined with the phonopy package (32) on supercells of 288 
(MgH 2 ), 160 (MgH 4 ), 351 (MgHi 2 ) and 272 (MgHi 6 ) atoms. 
The electronic densities of states and phonon band structures 
obtained with VASP for MgH 2 , MgH 4 and MgHi 2 at 180, 
100 and 140 GPa showed good agreement to those computed 
using the Quantum Espresso (QE) program [ 33 ], and the com- 
putational settings described below. 

In the QE calculations the H and Mg pseudopotentials, ob- 
tained from the QE pseudopotential library, were generated 
by the method of von Barth and Car with Is 1 and 3s 2 3p° 
valence configurations, along with the Perdew-Zunger local 
density approximation [34]. We choose this particular Mg 
pseudopotential as it has been employed in numerous lattice 
dynamical studies of materials under pressure, see for exam- 



ple Ref. (35). Plane wave basis set cutoff energies, which 
gave an energy convergence to better than 0.3 mRy/atom were 
55, 75 and 90 Ry for MgH 2 , MgH 4 and MgH i2 , respec- 
tively, and al6 x 16 x 16 Brillouin-zone sampling scheme 
of Methfessel-Paxton (36) with a smearing factor of 0.02 Ry 
was employed. Density functional perturbation theory, which 
is implemented in QE, was used for the phonon calculations. 
The electron-phonon coupling matrix elements were calcu- 
lated using al6xl6xl6 /c-mesh and 4x4x4 g-mesh 
for MgH 2 and MgH4, along with al2 x 12 x 12 /c-mesh and 
3x3x3 g-mesh for MgH i2 . The electron phonon coupling 
(EPC) parameter, A, was calculated using a set of Gaussian 
broadenings in steps of 0.005 Ry. The broadening for which A 
was converged to within 0.005 was 0.035, 0.025 and 0.040 Ry 
for MgH 2 , MgH 4 and MgH i2 , respectively. The supercon- 
ducting transition temperature, T c , has been estimated using 
the Allen-Dynes modified McMillan equation [ 37 ] as: 
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where uj\ og is the logarithmic average frequency and /i* is the 
Coulomb pseudopotential, often assumed to be between —0.1- 
0.13. 

The molecular calculations on H^ - and MgH i2 were per- 
formed using the ADF software package (38) and the revPBE 
gradient density functional. The band structures of select 
phases (see the Supporting Information, SI) were calculated 
using the tight-binding linear muffin-tin orbital (TB-LMTO) 
method (39), the VWN local exchange correlation potential 
along with the Perdew-Wang GGA. 



III. RESULTS AND DISCUSSION 
A. Squeezing MgH 2 

The potential for reversible hydrogen storage has resulted 
in much interest not only in crystalline MgH 2 , but also in 
nanoparticles based on the structure of the various phases 
of this solid (40H44) . At 1 atm a-MgH 2 assumes a Ti0 2 - 
rutile-type geometry and it undergoes a series of structural 
transitions, a (P^/mnm) — » 7 (Pbcn) — » f3 (Pa3) — » 
S (Pbc2 1 ) e (Prima), at 0.39, 3.84, 6.73 and 10.26 GPa 
(45l|46). At around 165 GPa the Prima phase slightly distorts 
to a higher symmetry P^jmmc Ni 2 In-type configuration, 
£-MgH 2 , see Fig. [I] which was shown to be dynamically sta- 
ble at 180 GPa. At lower pressures BaH 2 (43118), CaH 2 (49), 
and SrH 2 [ 50 ] also adopt this structure. However, whereas the 
heavier alkaline earth dihydrides with P6 3 /mmc symmetry 
are insulating, in MgH 2 metallization is predicted to occur by 
170 GPa at the PBE-level of theory (5T) . Our evolutionary 
runs did not reveal any other phases of MgH 2 up to 300 GPa, 
a pressure at which the metallicity in (-MgH 2 persists. Calcu- 
lations using the HSE06 screened hybrid functional [ 27 ] con- 
firm the metallicity at -180 GPa (see the SI). 

£-MgH 2 is metallic as a result of the closure of an indi- 
rect band gap. A flat band displaying H~ s-character rises 
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above the Fermi level, Ep, around T — A, and a steep band 
that boasts primarily Mg s and a little bit of Mg p-character 
falls below E F around the i^-point, as illustrated in Fig. [T] 
(-MgH2 contains half the number of valence electrons per for- 
mula unit as does MgB 2 , which becomes superconducting at 
39 K [52]. The B a- and H s-bands are comparable, in partic- 
ular the "holes" at the top of the band spanning from T to A. 
However, whereas MgB2 consists of hexagonal boron sheets 
with Mg 2+ intercalated in the hexagonal holes, the hexagonal 
network in (-MgH2 is made up of alternating Mg 2+ and H~ 
ions, with the second set of hydrides located in the hexagonal 
holes. So the Fermi surfaces arising from the B a- l53l and 
the H s-bands |51 ] are not identical. We also noticed some 
similarities between £-MgH2 and the most stable CsH phase 
above 160 GPa (54). Whereas the former can be thought of 
as layers of graphitic sheets of alternating Mg 2+ and H~ ions 
arranged in an ABABA... stacking with H~ sandwiched in 
between the two layers, in CsH half of the H~ have been re- 
moved so that only layers of Cs + and H~ are found. 

The density of states at the Fermi level, g(Ep), of £-MgH 2 
is 0.01 eV "Velectron at 170 GPa and by 300 GPa it decreases 
only slightly, see the SI. Despite the relatively low g(Ep) we 
were intrigued in the possibility of superconductivity in this 
system. Fig.[2]plots the phonon band structure and densities of 
states, phonon linewidths (7(0;)), Eliashberg spectral function 
(a 2 F(u)), and X(u) of C-MgH 2 at 180 GPa. 38% of the to- 
tal EPC parameter, A, is a result of the low-frequency modes 
below ~700 cm -1 which are associated primarily with the 
motions of the heavier Mg ions, whereas the region between 
750-2400 cm -1 , which is mostly due to motions of the H~ 
anions, contributes 62% towards the total A with the modes 
along T — K — M and H — L — A playing a dominant role. 
The modes above 2100 cm -1 are primarily caused by the set 
of H~ anions which are closer to the Mg 2+ cations (the Mg-H 
distances at 180 GPa measure 1.58 and 1.81 A). Despite the 
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FIG. 1: The electronic band structure of C~MgH2 at 180 GPa. The 
character of the bands which cross the Fermi level is highlighted. 
The Fermi energy is set to zero in all of the plots. A side and top 
view of ^-MgH2 is shown on the right, with the Mg/H atoms colored 
in red/ white. 
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FIG. 2: Phonon band structure, phonon density of states and the 
Eliashberg spectral function, a 2 F(uj), of C-MgH 2 at 180 GPa. Cir- 
cles indicate the phonon linewidth with a radius proportional to 
the strength. At this pressure A = 0.58, uJi og = 1111 K, and 
T c = 16 - 23 K assuming fi* = 0.13 - 0.1. 



modest g(Ep) the total EPC parameter is calculated as being 
0.58, and along with an uj\ og of 1111 K gives rise to a T c of 
16-23 K for /i* ranging from 0.13 — 0.1, respectively, via the 
Allen-Dynes modified McMillan equation. For comparison, 
the EPC parameter and uj\ og for BaH2 in the simple hexagonal 
structure at 60 GPa have been calculated as being 0.22 and 
780 K, respectively, giving rise to T c on the order of only a 
few mK R71 

Given the recent interest in compressed hydrogen-rich 
solids as potential superconductors, we began to wonder if H2 
may mix with MgH2 under pressure and what other metallic 
systems may be found. As will be shown in a moment, while 
the ( phase is preferred for MgH 2 above 165 GPa, it is not the 
most stable point on the hydrogen-rich Mg/H phase diagram 
at these pressures. 



B. Stabilization of the Polyhydrides 

The first ionization-potential of magnesium is ~2 eV larger 
than that of lithium, suggesting that stabilization of the mag- 
nesium polyhydrides will occur at a somewhat higher pressure 
than the lithium polyhydrides (that is, above 100 GPa fT3ll ). 
The ionic radius of Mg 2+ is about 30% smaller than that of 
Ca 2+ , and since CaH 4 was found to be the most stable stoi- 
chiometry at P = 50 - 150 GPa and CaH 6 at 200 GPa l20l . 
we may expect that the most favorable MgH n combination has 
2 < n < 6. How do these predictions — based upon chem- 
ical intuition that was developed by analyzing computational 
studies of compressed polyhydrides fT3jj20l — compare with 
the results from our evolutionary structure searches? 

The calculated enthalpies of formation, AHf, of the im- 
portant MgH n structures found in our searches are provided 
in Fig. [3] (see the SI for the plot of AHf vs. H2 composi- 
tion). At 100 GPa only MgH 4 is predicted to resist decom- 
position into MgH2 and H2. In fact, MgH 4 persists as hav- 
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ing the most negative AH F as the pressure is increased. By 
200 GPa all stoichiometries except for MgH2.5, MgH3 and 
MgH 3 .5 have a negative AH F . Because MgH 4 , MgHi 2 and 
MgHi6 lie on the convex hull at 200 GPa, which is provided 
in the SI, they are thermodynamically stable with respect to 
decomposition into other phases and MgH2/H2. Fig. [^illus- 
trates that whereas the AH F of MgH 4 starts to become neg- 
ative near 92 GPa, for MgHi 2 and MgHi 6 this occurs at 122 
and 117 GPa, respectively. Comparison of these findings with 
the predictions based upon our newly developed chemical in- 
tuition under pressure illustrates that we did reasonably well 
in predicting the stabilization pressures and most stable stoi- 
chiometry without performing any computations. Now, let us 
begin our exploration of the structural peculiarities and elec- 
tronic structures of the magnesium polyhydrides falling on the 
convex-hull. 



C. MgHzt: H 2 molecules and hydridic atoms 

The phase with the most negative enthalpy of formation 
throughout the pressure range studied, Cmcm-MgYi^ forms 
a distorted CsCl bcc lattice with H~ atoms taking up the ver- 
tex positions of the underlying cube-like cages. The body- 
centered positions are occupied by a 1:1 mix of Mg 2+ cations 
and H2 molecules with slightly elongated bonds, as high- 
lighted by the red and blue polyhedra in Fig. |4| a )- The size 
mismatch between the two results in the hydrides forming 
different sized cages around them. The same Cmcm-MgH 4 
structure was discovered in evolutionary runs carried out at 
100 and 200 GPa, and phonon calculations at 100 GPa showed 
it was dynamically stable. 

Interestingly, the lowest enthalpy 74/mmm-symmetry 
CaH 4 structure found using the particle-swarm optimization 
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FIG. 3: AH F of the MgH 4 , MgHi 2 , and MgHi 6 phases as a func- 
tion of pressure. The plot also illustrates the enthalpy of the various 
MgH2 phases with respect to the most stable structure at a given 
pressure, with the region below 15 GPa magnified in the inset. 
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FIG. 4: (a) A schematic illustration of the structural similar- 
ity between Cracra-MgFU (right) and the 74/rarara-symmetry 
CaH4 structure from Ref. |20| (left). The hydride cages around 
Mg 2+ /Ca 2+ and the H2 molecules are colored red and blue, re- 
spectively |56|. (b) The PBE electronic band structure of MgH4 at 
100 GPa. The character of the bands which cross the Fermi level 
is highlighted. Also shown on the right are the total electronic and 
hydrogen site-projected densities of states (DOS). 



technique 1 55 ] in a recent study l20l can be described the same 
way as Cmcm-MgH^. The difference between the two is the 
manner in which the cations and the H 2 molecules are dis- 
tributed: in CaH4 they are dispersed homogeneously, but in 
MgH4 they are arranged into interwoven zig-zag chains. Ad- 
ditionally, the hydride cages are more distorted in the MgH4 
structure with the smallest face which surrounds an Mg 2+ be- 
ing the one which links one magnesium cation to another. Ac- 
tually, the structure can also be viewed as sheets of MgF^ with 
H 2 molecules trapped in the larger gaps between the sheets, 
one such sheet is shown at the right of Fig.[4ja). 

Within PBE the total density of states of Cmcm-MgH 4 at 
100 GPa is quite small; in fact it is lower than that of the 
14/mmm analogue. The former is 42 meV/atom more stable 
than the latter at this pressure, but the volume of IA/mmm- 
MgH 4 is 0.04 A 3 per atom smaller than that of Cmcm- 
MgH4. The formation of the zig-zag like chains illustrated 
in Fig. |4ja) must therefore have a substantial impact on the 
electronic contribution to the enthalpy. 

As the pressure increases from 100 to 300 GPa, the short- 
est distance between the metal cations decreases from 2.56 A 
to 2.32 A, and the Mg 2p core bands broaden from 0.19 to 
0.58 eV due to core overlap. The shortest Mg-H - distance 
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decreases from 1.74 A to 1.55 A, as does the Mg-H 2 sep- 
aration (1.92 A/1.68 A). The molecular H2 bonds lengthen 
somewhat from 0.78 A to 0.79 A. This is slightly longer than 
the intramolecular distance in pure compressed H2 at these 
pressures, 0.73 A and 0.75 A, respectively. 

Since MgH4 contains H~ units, it becomes metallic as a re- 
sult of pressure induced band broadening and eventual over- 
lap, like LiH2 and NaHg. Metallization occurs within PBE 
already by 20 GPa, even though Cmcm-MgH 4 is not stable 
with respect to decomposition into H2 and MgH2 at this pres- 
sure. A band displaying H~ s-character rises above ^ at 
the F-point, whereas bands exhibiting predominantly Mg sp 
and a little bit of H2 a* -character fall below it at the Z and 
T-points (see Fig.|4jb), and the "fat bands" in the SI). Calcu- 
lations using the HSE06 functional increase the metallization 
pressure to ~150 GPa (see the SI). 

D. MgHis and MgHi 6 : the belt 

The absence of H~ anions in the other structures lying on 
the convex hull at 200 GPa, MgHi2 and MgHi6, has impor- 
tant consequences for their electronic structure. The hydro- 
genic sublattices of these phases are composed solely of H^ - 
molecules with slightly stretched bonds, just like the CaHi 2 
structure predicted in a previous study l20l . The presence of 
hydridic hydrogens in systems with a small mole fraction of 
hydrogen, and their absence in phases with a high hydrogen 
content was rationalized by Wang and co-workers by con- 
sidering the formal number of 'effectively added electrons' 
(EAE) which are donated from the alkaline earth metal va- 
lence bands into the H 2 a*-bands [20] for various CaH n stoi- 
chiometrics. If the EAE is small, the dihydrogen bond simply 
stretches as a result of the population of H 2 cr*, and hydridic 
hydrogens are not formed. If the EAE is large enough, the 
molecule dissociates into two H~ units. 

The most stable MgHi2 structure found in our evolution- 
ary searches had i?3-symmetry and was shown to be dy- 
namically stable at 140 GPa. The hexagonal building block 
for this phase consists of an Mg 2+ cation surrounded by six 
H^" 1 ^ 3 molecules, as illustrated in Fig. [5Ja). The dihydro- 
gen molecules in these hexagons appear to arrange in a 'belt' 
around the metal cation in a side-on fashion. The hexagons 
are tiled in parallel sheets forming an ABC ABC... close- 
packed structure. This description is, however, somewhat mis- 
leading as the distance between Mg 2+ and a hydrogen atom 
within the MgHi2 building block is comparable to the short- 
est distance between the metal ion and a hydrogen atom in a 
different building block. For example, the intra- and inter- 
building block Mg-H measurements are 1.76 and 1.81 A, re- 
spectively, at 140 GPa. The dihydrogen bonds in MgHi2 are 
substantially elongated to 0.82 and 0.81 A at 140 and 300 GPa, 
when compared with the typical bond lengths of pure H2 at 
these pressures (0.74/0.75 A). A gas phase geometry opti- 

— 1/3 

mization on a H 2 / molecule at 1 atm showed that the bond 
stretched from 0.75 to 0.79 A. So the expansion of the H-H 
bond observed in MgHi 2 is consistent with the partial filling 
of the H 2 a*-bands. 




FIG. 5: (a) A 2 x 2 x 2 supercell of i?3-MgHi 2 where the polyhedral 
units are highlighted to emphasize the ABC ABC... close-packed ar- 
rangement. The MgHi2 building blocks are colored such that the red 
polyhedra all lie in the same plane, (b) A supercell of P — l-MgHi6. 
The MgHi2 building blocks are colored to emphasize the similarity 
to the structure in (a). The blocks colored red or green lie in the same 
plane, and the hydrogens colored in blue do not belong to an MgHi2 
building block |56l . 



The lowest enthalpy MgHie configuration we found, which 
exhibited P — 1 symmetry and was dynamically stable at 
130 GPa, could be thought of as being made up of similar 
MgHi2 units with excess H2 molecules stuffed between the 
hexagons, as illustrated in Fig.j5jb). These 'interstitial' dihy- 
drogens lie on top of the Mg 2+ cations in the MgHi 2 building 
blocks. The intramolecular distance in H^ 4 at 1 atm was 
calculated as being 0.78 A, whereas the bond lengths of the 
dihydrogens in MgHi 6 range from 0.76 to 0.85 A at 130 GPa, 
with the interstitial hydrogen molecules being shorter (0.76- 
0.78 A) than those in the belt (0.79-0.85 A). The unequal bond 
lengths in the H2 units may be a result of different charge 
states, or be due to their local bonding environments, or both 
1571 . As the pressure increases, the intramolecular H-H bond 
lengths decrease, with the interstitial and belt hydrogens mea- 
suring 0.74-0.76 A and 0.78-0.84 A, respectively, at 300 GPa. 

Since MgHi2 and MgHi6 are metallic as a result of the par- 
tial filling of the H 2 a*-bands like LiH 6 , they have a high 
density of states at Ep, see Fig. [6] Both phases remain good 
metals upon compression up to at least 300 GPa, and the Mg 
2p bands broaden only slightly to ~0. 15-0.2 eV at 300 GPa as 
a result of core overlap. A comparison of the DOS of MgHi 2 
at 140 GPa calculated with PBE and the HSE06 screened hy- 
brid functional showed that latter valence DOS was slightly 
broader (as expected for metallic systems (58l [59)), and the 
core Mg 2p bands shifted to lower energies. However, the 
g(Ep) computed with the two functionals was essentially 
the same (see the SI). This is in-line with our previous re- 
sults which showed that g(Ep) is relatively insensitive to the 
choice of the functional in sodium polyhydrides that did not 
contain hydridic hydrogens lfT4l . The nearly free-electron like 
valence bands of MgHi 2 and MgH i6 exhibit primarily H s- 
character, with a little bit of Mg sp spanning throughout (see 
the fat bands in the SI). Nonetheless, the computed DOS of 
these structures match quite well with a hypothetical system 



> 

CD 



0.12 



0.1 



0.08 



0.06 



0.04 



0.02 



-30 




(a) 



-20 -10 
Energy (eV) 



20 -10 
Energy (eV) 



FIG. 6: The PBE valence densities of states of (a) MgHi2, and (b) 
MgHi6. Note that the Fermi level lies just below (directly after) a 
sharp peak in the DOS in MgHi 2 at 140 GPa (MgHi 6 at 130 GPa). 



where the metal cations have been removed, the (Hi 2 ) 2_ and 
(Hi6) 2_ lattices (see the SI), suggestive of almost full ioniza- 
tion of the Mg valence 3s electrons into the H 2 <r*-bands. 



E. The MgHi2 Building Blocks 

A number of MH12 clusters, where M is a transition metal 
atom, have been predicted as being stable in the gas phase by 
quantum chemical calculations [60]. Moreover WH 4 (H 2 )4, 
which contains four hydridic hydrogens and four dihydrogens 
bonded to the metal center in a side-on fashion, has been 
made in a neon matrix 16T1 . Unsurprisingly, molecular calcu- 
lations on both the optimized MgHi 2 building blocks as well 
as (MgHi 2 ) 2+ with Dq^ symmetry revealed numerous imagi- 
nary frequencies, so these clusters will not be stable at 1 atm. 
Nonetheless, it is instructive to compare the electronic struc- 
ture of the MgHi2 building block with the DOS computed for 
J R3-MgHi2 at various pressures. 

Whereas in an optimized MgHi2 cluster the Mg-H and 
intramolecular H-H distances were calculated as being 2.20 
and 0.79 A, these values measure 2.34-2.58/1.76-1.78 A and 
0.79/0.82 A in R3-MgU 12 at 5/140 GPa. So the gas phase 
cluster is actually a little bit more 'compressed' than the build- 
ing block within the solid at 5 GPa. The energy level diagrams 
and canonical molecular orbitals (MOs) of these clusters are 
illustrated in Fig.|7Ja). The six MOs lowest in energy resem- 
ble the canonical MOs of benzene, except that they do not 
contain a node which bisects the dihydrogen molecules. A 
large gap is found between the highest occupied molecular 
orbital (HOMO) and HOMO-1 of the optimized cluster, but a 
small gap separates the HOMO and lowest unoccupied molec- 
ular orbital (LUMO). The HOMO-HOMO-1 gap decreases as 
the cluster is compressed, as do the analogous sets of bands in 
the extended system, as illustrated in Fig.|7Jb). 
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FIG. 7: (a) Calculated molecular orbital level diagram of an MgHi2 
cluster when optimized in the gas phase at 1 atm (left) and in the 
geometry which it adopts in the i?3-MgHi2 solid at 5 GPa (center) 
and 140 GPa (right). The energies of the MOs at 140 GPa are given 
to the right in eV. (b) Calculated site-projected densities of states of 
#3-MgHi2 at 5 and 140 GPa. 



The frontier orbitals of the 1 atm cluster contain substantial 
character arising from H, and the HOMO also has an impor- 
tant contribution from the Mg 3p and the LUMO from the Mg 
3s orbitals. Our fragment orbital analysis [62] shows that the 
LUMO+1 contains primarily H-character with less than 4% 
of Mg 3d states mixed in. In the 'less-compressed' cluster 
extracted from the solid at 5 GPa the frontier orbitals swap 
positions so that the HOMO displays Mg 3s and the LUMO 
Mg 3p-character. This is in-line with the projected densi- 
ties of states which illustrate that the predominant contribution 
around the Fermi level is due to Mg s at 5 GPa. At 140 GPa the 
gap between the hydrogenic and the metallic bands between -6 
to -3 eV closes in the solid, and the HOMO and HOMO-1 or- 
bitals in the cluster come closer together in energy. Moreover, 
as a result of the pressure induced s — » p transition in Mg, the 
states around the Fermi level contain about an equal amount 
of Mg s and p-character. In the cluster the energy ordering 
of the molecular orbitals also changes, so that the HOMO is 
more Mg p-like, and the LUMO displays primarily H s char- 
acter. So the gas phase clusters are able to mimic some of the 
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FIG. 8: Phonon band structure, phonon density of states and the 
Eliashberg spectral function, a 2 F(uo), of MgH4 at 100 GPa. Cir- 
cles indicate the phonon linewidth with a radius proportional to 
the strength. At this pressure A = 0.74, oji og = 941 K, and 
T c = 29 - 37 K assuming /x* = 0.13 - 0.1. 



FIG. 9: Phonon band structure, phonon density of states and the 
Eliashberg spectral function, a 2 F(u), of MgHi 2 at 140 GPa. Cir- 
cles indicate the phonon linewidth with a radius proportional to 
the strength. At this pressure A = 0.73, cji og — 1554 K, and 
T c = 47 - 60 K assuming /x* = 0.13 - 0.1. 



essential features of the projected DOS of the solid at different 
pressures. 



F. Superconductivity in the polyhydrides 

Our computations have estimated the T c of the classic al- 
kaline earth hydride, P6s/rarac-MgH 2 , as being ~20 K at 
180 GPa. We wondered what the T c of the aforementioned 
hydrogen-rich phases would be in comparison, and if T c 
would be influenced by the hydrogenic sublattice of the poly- 
hydride. Computations have been carried out on Cmcm- 
MgH 4 and R3-MgU 12 at 100 and 140 GPa, respectively, 
since these are pressures slightly larger than when AHp was 
computed as becoming negative, and since both systems were 
found to be metallic within PBE at these pressures. 

Within PBE the density of states at the Fermi level of 
Cmcm-MgH 4 at 100 GPa is comparable to the one calcu- 
lated for C-MgH 2 at 180 GPa. The total EPC parameter is 
nearly 28% larger, however. In both phases the low-frequency 
modes which are mostly due to the vibrations of the Mg atoms 
contribute ~0.2 to the total A; compare Fig.[8]and Fig. [2] The 
main reason why the overall EPC parameter is larger for the 
polyhydride than the classic hydride is a result of the total 
coupling provided by the modes between ^500-2500 cm -1 , 
which are primarily due to the motions of the hydrogen atoms. 
Whereas the classic hydride contains only hydridic hydro- 
gens, there is an equal number of H2 and H~ hydrogens in 
MgH 4 . The H 2 vibron, located around 3500 cm -1 contributes 
only 1.5% to the total coupling strength in the polyhydride. 
This is not surprising, since the bands crossing the Fermi level 
displayed only a small amount of H2 s-character. Despite the 
higher A that MgH 4 has as compared to MgH 2 , the prefactor 
in the modified McMillan equation, u\ og , is calculated as be- 
ing 15% smaller so the total T c of MgH 4 is estimated as being 
~14 K higher than that of the classic alkaline earth hydride. 



Despite the increased hydrogen content and higher density 
of states at the Fermi level, the total EPC parameter of R3- 
MgHi2 at 140 GPa is almost the same as of Cmcm-MgH 4 
at 100 GPa. However, a comparison of the X(uj) in Fig. [9] 
and Fig. [8] shows that the relative contributions to the overall 
A are quite different in the two polyhydrides. MgHi2 does 
not contain any hydridic hydrogens, and is metallic because 
of the partial filling of the H2 a* -bands. In contrast to what 
was found for MgH 4 , the H 2 vibron contributes about 11.1% 
to the total A. This corresponds quite well to the 11.3% cal- 
culated for a compressed KH 6 phase [|T9l whose hydrogenic 
sublattice only contained H^ - motifs. The EPC associated 
with the low-frequency modes below 400 cm -1 , which are 
dominated by the motions of the heavier metal atoms, is about 
a third of the amount calculated for the phases containing H~ 
units. The main contribution to A, 79%, arises from the inter- 
mediate frequency regime, which is primarily due to the H 2 
motions. The reason why the T c of MgHi 2 is estimated as be- 
ing ~20 K higher than that of MgH 4 is due to the larger cji 0g . 
Unfortunately, the computational expense precluded us from 
calculating the EPC parameter of compressed P — l-MgHi6, 
or from exploring the pressure dependence of T c . 

The values we calculate for the total EPC parameter and 
the T c of the magnesium polyhydrides falls within the range 
of 0.5-1.6 and 10-139 K, respectively, computed for a number 
of hydrogen-rich systems l20l . We show that comparable A 
values may be obtained for polyhydrides with very different 
hydrogenic sublattices, but their u\ og and therefore T c may 
differ. The magnesium polyhydrides are predicted to have a 
larger T c than MgH 2 under pressure, and phases with a larger 
mole percent ratio will likely have a higher T c . 
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IV. CONCLUSIONS 

Evolutionary structure searches coupled with density func- 
tional theory calculations are used to predict the most sta- 
ble structures and stoichiometries of the magnesium polyhy- 
drides, MgH n with n > 2, underpressure. The thermodynam- 
ically stable structures found in this study have a hydrogenic 
sublattice containing H~ anions and H 2 units (MgH 4 ), or H 2 
molecules which are less strongly bonded than those found in 
pure molecular hydrogen at 1 atm (MgHi 2 and MgHi 6 ). 

Metallization in P63 / rarac-MgH2 occurs as a result of an 
H~ s-band rising above, and a Mg sp band falling below 
the Fermi level. T c is estimated as being between 16-23 K 
at 180 GPa, with a sizable contribution to the total electron 
phonon coupling parameter arising from vibrations related to 
both the hydrogen and magnesium atoms. 

MgH4, which starts to become thermodynamically sta- 
ble with respect to decomposition into MgH 2 and H 2 near 
100 GPa is found to contain one H 2 molecule and two hydridic 
hydrogens per Mg 2+ cation. Metallization occurs as a result 
of pressure-induced band gap closure, but the density of states 
at the Fermi level is quite low. Around 120 GPa other stoi- 
chiometries, whose hydrogenic sublattices contain only H^ - 



molecules with slightly stretched bonds, emerge as being ther- 
modynamically stable. MgHi 2 and MgHi6 are metallic in 
part as a result of the partial filling of the H 2 a* -bands and 
have a high density of states at the Fermi level. Their elec- 
tronic structure at various pressures can be traced back to the 
molecular orbital diagram of their building block, the MgH i2 
cluster. Despite the very different hydrogenic sublattices, both 
MgH4 and MgHi 2 are found to have similar electron phonon 
coupling parameters. The main reason why the T c of MgH i2 
at 140 GPa is calculated as being larger than that of MgH4 at 
100 GPa, 47-60 K vs. 29-37 K assuming typical values of /A 
is because of the larger average logarithmic frequency com- 
puted for MgH i2 . 
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